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Abstract. The spectrum of (3 CMa (Mirzam) between 1000 and 1200A was recorded 

I Profile Spectrograph (IMAPS) during its orbital flight on the ORFEUS-SPAS II mission 
\ in 1996. New information about interstellar absorption lines of Cn, Cn*, Ni and Ol 
from this spectrum are combined with the HST results reported by Dupin & Gry (1998) 
^ ■ to arrive at new conclusions about the physical properties of the absorbing gas in front 
' of j3 CMa. For two prominent velocity components centered at heliocentric velocities 

H ; 

of +20.0 and +30.5 km s , designated by Dupin & Gry as Components C and D, re- 
spectively, we use determinations of iV(Cn*)/iV(Cll) and iV(Mgl)/iV(MgIl) to solve for 
temperatures and electron densities. From our knowledge that oxygen and nitrogen have 
their ionizations coupled to that of hydrogen through charge exchange reactions, we can 
derive the hydrogen ionizations by comparing these elements to sulfur, which is likely not 
to be depleted onto dust grains. For Component C with an approximate column density 
of neutral and ionized hydrogen iV(Htotai) = 6 x 10 18 cm -2 , we find that the neutral frac- 
tion n(Hi)/n(H to tai) = 0.25, 400 < T < 6500 K, and 0.08 < n(e) < 0.6 cm- 3 , while for 
Component D with iV(H to tai) = 1-2 x 10 19 cm -2 , we arrive at n(Hi)/n(H to tai) = 0.035, 
8000 < T < 14 000 K, and 0.09 < n(e) < 0.2cm- 3 . The relatively large ionization frac- 
tions of H can arise if the clouds are about 130 pc away from us, so that they are exposed 
to the strong, ionizing radiation fields from e and (3 CMa. The presence of Si in indicates 
the existence of additional gas with even higher levels of ionization. 
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1. Introduction 

The line of sight toward the Bl II-III star CMa in the direction (In, 6//)=(226?, -14?) 
has raised interest since its first study in the UV from Copernicus observations when Gry 
et al. (1985) showed that it contained very little neutral gas and that a great majority 
of the material on the line of sight was ionized. Indeed for a distance now set to 153 pc 
by Hipparcos measurements, Gry et al (1985) determined that iV(Hi) was somewhere 
within the range 1 — 2.2 10 18 cm -2 based on the interstellar Ly/3 absorption profile in 
the spectrum of (3 CMa recorded by Copernicus, and showed that the total (neutral 
and ionized) hydrogen column density was ten times higher using the combined abun- 
dances of Sn and S in compared to that of Hi multiplied by the cosmic ratio of S to H. 
EUVE observations independently confirmed both findings. Cassinelli et al (1996) found 
N(Hl)=2.0 ± 0.2 10 18 cm~ 2 by fitting the Lyman limit absorption in the EUV spectrum 
of the star. Evidence for ionization comes from the lower limits for the continuum ab- 
sorption by He I in the spectrum of (3 CMa which indicate that iV(Hei) > 6. 10 17 cm~ 2 
(Aufdenberg, et al. 1999) or > 1.4 10 18 cm~ 2 (Cassinelli, et al. 1996), depending on as- 
sumptions about the drop in stellar flux across the Hei ionization edge. The fact that 
iV(Hel) 3> O.liV(Hi) indicates that a substantial fraction of the hydrogen toward (3 CMa 
is ionized. 

The weak presence of H I gas in this region of the sky is particularly evident in the 
neighbouring line of sight toward e CMa (In — 239?8, bn = — 11?3) where the Lyman 
limit absorption is even lower (Cassinelli et al, 1995). Gry et al (1995) found a Hi column 
density upper limit of 5 10 17 and showed furthermore that almost all the neutral gas is 
very close to the Sun, distributed in two small components which are also detected 
with similar column densities in the line of sight to Sirius 2.7 pc away (Lallement et 
al.,1994). Because of the scarcity of Hi in this region, e CMa and (3 CMa are the two 
strongest sources of EUV radiation, and the two stars, especially e CMa, dominate the H- 
ionization field close to the Sun (Vallerga & Welsh 1995 ; Vallerga 1997). The conditions 
along these lines of sight, especially their ionization structure, strongly influence the 
nature and ionization of the local interstellar matter. Yet even with this in mind, we find 
that the ionization of the local ISM is not very well understood. For instance, with the 
radiation from e and (3 CMa plus other easily identifiable sources of ionizing radiation, it 
is difficult to reconcile the fact that helium appears to be more ionized than hydrogen, 
as revealed by evidence from the EUVE spectra of several white dwarfs (Dupuis et al. 
1995). 
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High resolution observations with GHRS on board HST allowed Dupin and Gry (1998) 
to show that the bulk of the gas in the (3 CMa sight-line is distributed in two main 
components separated by 10 km/s, which are warm, only slightly depleted and both 
mostly ionized. However the analysis of the ionization processes was hampered by the 
fact that several species and especially the neutral species had been observed with a 
limited resolution of about 20 000 which did not allow for a clear separation of the two 
components. The relative distribution of matter in the two components therefore had 
uncertainties that allowed only lower limits to the ionization fractions to be derived. 

Here we present new observations performed with the Interstellar Medium Absorption 
Profile Spectrograph (IMAPS) (Jenkins, et al. 1996), with the aim of gaining further 
insights on the ionization structure and various physical processes in the local interstellar 
medium. The higher resolution of 60 000 provided by IMAPS allowed us to resolve the 
absorptions of the two main components which in turn permitted the derivations of 
ionization fractions in each cloud. We also determined the clouds' electron densities using 
the observation of Cn*, which ultimately led to new insights on the ionization process. 

2. Observations and data reduction 

The observations were carried out by IMAPS when it was operated on the ORFEUS- 
SPAS II mission that flew in late 1996 (Hurwitz, et al. 1998). IMAPS is an objective- 
grating echelle spectrograph that was designed to record the spectra of bright, early-type 
stars over the wavelengths from ~ 950 A to ~ 1150 A with a high spectral resolution. 
For more details on the instrument see Jenkins et al. (1998, 1996). 

The spectra were extracted from the echelle spectral images using special procedures 
developed by one of us (EBJ) and his collaborators on the IMAPS investigation team. 
A zero point in the wavelength calibration was made by measuring the O I* and O I** 
telluric lines detected at 1040.943 A and 1041.688 A. Table |] presents the interstellar 
lines detected in the data, with wavelengths and /-values taken from Morton (1991). 

3. Analysis 

The structure of velocity components in the sight-line toward /3 CMa has already been 
determined by the analysis of GHRS HST data (Dupin & Gry 1998) : four components 
were detected, two of which (designated as Components C and D) were strongly dom- 
inant, and the remaining two (A and B) showed up only for the strongest lines. In the 
study conducted by Dupin & Gry (1998), only the lines between 1800 A and 3000 A 
(Fen, Mgn, Mgi and Sin lines) were observed at a resolution A/AA ~ 80 000. Because 
of the failure of the Echelle A grating at the time of the observations, all the species that 
had lines between 1150 A and 1800 A were observed at medium resolution (~ 20 000). 
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Table 1. Interstellar absorption lines detected in the (3 CMa IMAPS data. 



Element 


A 




f 
I 




Cn 


1036.337 


1.23 


10" 


- 1 


Cn* 


1037.018 


1.23 


10" 


- 1 


Ni 


1134.165 


1.34 


10" 


-2 




1134.415 


2.68 


10" 


-2 




1134.980 


4.02 


10" 


-2 


Nil 


1083.990 


1.03 


10" 


-1 


Oi 


1039.230 


9.20 


10" 


■3 


Sin 


1190.416 


2.50 


10" 


-1 




1193.290 


4.99 


10" 


-1 



Table 2. Velocities and column densities for various species toward (3 CMa. 



Component v e iV(Cn) JV(Cn*) AT(Ni) N(Oi) 

(kms- 1 ) (cnr 2 ) (cm" 2 ) (cm" 2 ) (cm" 2 ) 

A 3.0 ~ 1.7 10 13 < 1. 10 11 < 1. 10 11 < 1. 10 11 

B 11.0 ~2 10 15 < 1. 10 11 <4. 10 12 7.5 ±2.5 10 

C 20.0 saturated 2.4±0.110 13 9.6 ± 0.4 10 13 7.6 ± 0.4 10 

D 30.5 saturated 1.2±0.2 10 13 2.7±0.3 10 13 2.2 ± 0.2 10 



With a resolution of ~ 60 000, the IMAPS spectrograph now allows us to analyse lines of 
some important species like N I and O I at a resolution comparable to that of the GHRS 
Echelle data. As a consequence, we can resolve partially clouds C and D which are sep- 
arated by only 10 kms- 1 . To facilitate a comparison between the IMAPS data and the 
HST data that were presented by Dupin & Gry (1998), we interpreted the column density 
vs. velocity profiles in terms of the components defined by them. However, the data are 
not precise enough to rule out other, slightly different combinations of components that 
could produce fits that are just as acceptable. 

The relative velocities between the four clouds are precisely known from the Dupin 
& Gry (1998) study. The heliocentric radial velocities found in our study are larger than 
those of Dupin & Gry (1998) by about 5.5 kms" 1 . This difference is larger than the stated 
error for the GHRS velocity zero points, but we favor the velocities from the IMAPS data 
because we could use the telluric absorption lines as a reference. The velocity of C n* is 
especially good since it is very close to the O I** line at 1041. 688A in the spectral format 
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of the IMAPS echelle and cross disperser gratings. The velocity dispersions b for C II, N I 
and O I in Components C and D are not well determined, since some blurring of the lines 
is caused by the instrument. However, they should be constrained by fe(Mgii) on the low 
end and the measured widths on the upper end. Thus, the principal free parameters are 
the column densities of various species in cloud C and cloud D. For the strong lines of N I 
and O I we see some evidence for the presence of Component B, and we can place upper 
limits for the strength of A. 

The principal source of uncertainty for the column densities is the determination of 
the background intensity caused by scattered light from the echelle and cross disperser 
gratings in IMAPS. As is evident from Jenkins, et al. (1996 - see their Fig. 12) there is 
some overlap of energy from adjacent echelle orders, precluding the use of the interorder 
regions for measuring the background levels. We could determine the overall behavior of 
the scattered light by examining the trends for the bottoms of strong stellar absorption 
features. For the background level below the line of C n* at 1037. 018A, we could use the 
nearby, saturated absorption from interstellar Cn at 1036. 337A. 

Fig. Q shows the features recorded by IMAPS after they were normalized to the stellar 
continuum. Overplotted on these data are the predicted absorption profiles for the column 
densities and central velocities given in Table |2| Selected results for Components C and 
D derived by Dupin & Gry (1998) are listed in Table ||, as an aid to following the 
interpretations that we present in §|| We made no attempt to measure the Si n lines 
listed in Table |l|, since they were saturated and not nearly as useful as the transition at 
1808 A observed by Dupin & Gry (1998). Conversely, our Nl lines are not saturated, but 
the weakest member of the 1200 A multiplet observed by Dupin & Gry had a central 
optical depth To = 1.7 for Component C [derived from our jV(Ni) and their b value]. The 
lack of saturation and higher resolution for the IMAPS recordings of N I absorption lead 
to our preferring the values derived here to those of Dupin & Gry, even though there is 
some uncertainty in the IMAPS background level. 

The O I transition at 1302 A observed by Dupin & Gry was hopelessly saturated, 
which led to their only being able to express a lower limit for N(Oi). Fig. [l] shows that 
the IMAPS recording of the much weaker transition at 1039 A shows a fairly strong 
saturation. If there are narrow, unresolved subcomponents within Component C that 
are far more saturated than what we see in the apparent residual intensity of the line, 
our listed column density may be below the true value. In principle, one could sense the 
presence of such components by observing that the N I absorptions do not grow as fast as 
their increases in transition strengths. However the lines that we observed here are too 
weak to show this effect well. Nevertheless, it is interesting to note that by observing the 
strong multiplet of Nl at 1200 A Dupin & Gry (1998) found AT(Ni) = 8.1 ±0.8 10 13 cm~ 2 



6 



Jenkins, Gry & Dupin: The Sight-line toward f3 CMa 




1134.15 1134.20 1134.25 1134.30 1134.35 1134.40 
Wavelength 



I I I 1 I I I I 1 

¥rf\M» 1 


i i i i i i i i i i i i 


: o i W 





039.20 1039.25 1039.30 1039.35 1039.40 
Wavelength 




1134.40 1134.45 1134.50 1134.55 1134.60 
Wavelength 




1036.35 1036.40 1036.45 1036.50 
Wavelength 



1 1 1 1 1 1 1 1 1 1 1 1 1 


i i i i i i i i i i i i i 


; N I 

i t t t. t i i i i i ] i i 





0.0 

1134.95 1135.00 1135.05 1135.10 1135.15 1135.20 
Wavelength 

1.2 r 




1037.00 1037.05 1037.10 1037.15 1037.20 
Wavelength 



Fig. 1. Interstellar lines observed with IMAPS. The histogram-style tracings represent 
the recorded data after normalization to the stellar continuum, and the smooth lines 
show the fits obtained from our analysis. The tick marks above the absorption features 
denote (left to right) the positions of the velocity Components B, C and D identified in 
the text, except for Cn which shows Component A (far left) in addition to the others. 



for Component C (including the hidden local clouds )[], which is only 0.84 times the value 
that we obtained (see Table |2|). (Our and their results for Component D agree however.) 
Thus, when lines reach a strength similar to that of the O I line at 1039 A, the column 
densities might be understated by a factor of about 0.84. 



4. Interpretation 

1 Note that in Dupin & Gry Table 2, the listed column densities for Component C had a 
contribution from the LIC and the other local cloud detected in both a CMa and e CMa 
sightlines subtracted off. This contribution is negligible for all species but O I and N I for which 
it represents about 25% of the total absorption in Component C. 
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Table 3. Selected results from Dupin & Gry's (1998) study of (3 CMa. 
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Species Comp. C Comp. D 

iV(Sll) (cm" 2 ) 1.1±0.4 10 14 2.3±0.5 10 14 

AT (Si n) 5.8 ± 0.5 10 13 7.5 ± 0.3 10 13 

TV(Siin) <5 10 12 (1.5-10)10 14 

iV(Mgi) 8.0 ± 0.6 10 10 2.5 ± 0.3 10 11 

iV(Mgn) 3.1 ± 0.6 10 13 5.9 ± 0.8 10 13 

iV(Feii) 1.2±0.110 13 1.0±0.110 13 



4--1- Radial velocities 

As noted earlier, there is a systematic offset between the velocities reported here and those 
given by Dupin & Gry (1988). Our new velocity for Component C at 20kms _1 makes 
it consistent with that of gas in the Local Interstellar Cloud (LIC), moving at a velocity 
of 26kms~ 1 towards the direction I = 186° and b = —16° (Lallement & Bertin 1992), 
whereas Dupin & Gry felt that the LIC was buried in the gap between Components C and 
D. However, there are two important problems for identifying the origin of Component C 
with the LIC. First, the column densities of Mg n and Fe II in Component C (see Table ||) 
are considerably larger than A(Mgn) = 1.65 10 12 cm -2 and TV(Fen) = 8.5 10 11 cm -2 
found in the LIC component at 18.8 kms -1 in the spectrum of a CMa (Lallement, et 
al. 1994) which is only 5.5° away in the sky, and than TV(Mgn) = 3.0 10 12 cm~ 2 and 
iV(Feii) = 1.35 10 12 cm~ 2 found in the LIC component at lTkms" 1 in the spectrum of 
e CMa (Gry, et al. 1995), 17° away in the sky. Also, iV(Hl) = 5 10 17 cm- 2 in front of 
a CMa (Holberg, et al. 1998) is only about a third as large as the Hi column density that 
we identify with Component C [iV(Hl) = 1.6 10 18 cm~ 2 , taken from iV(Ni) and A(Oi) 
after correcting for depletion - see § |4.2| below.] Even though a CMa is only 2.6 pc away, 
it seems to lie beyond the boundary of the LIC because Lallement, et al. (1994) detected 



another (somewhat weaker) component at v = lSkms^ 1 . Second, we show in §4.5 that 
the ionizing radiation field that is needed to sustain the ionization is considerably higher 
than that found for the local vicinity. Thus, the agreement of velocity of Component C 
with that of the LIC is probably just coincidental. 
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4-2. Column densities, depletions, and ionization fractions 

We define a depletion factor 5(X) of an element X in terms of its column density N(X), 
relative to that of hydrogen A^H), by the expression 



where (H/JT) cosm i C is the cosmic abundance ratio.g For such a depletion to be meaningful 
in any given circumstance, one must account for differences in the observed ionization 
stage(s) of element X relative to that adopted for H (either Hi or Htotai)- We address 
this consideration in the cases that follow, starting with the simplest ones involving the 
elements nitrogen, oxygen and sulphur. 

The relative ionizations of N I and O I are closely coupled to that of H I through 
the large charge exchange reaction rates that arise from their nearly equal ionization 
potentials (Field & Steigman 1971; Butler & Dalgarno 1979). In the diffuse interstellar 
medium, these two elements are only mildly depleted below their cosmic abundances 
relative to hydrogen (Hibbert, et al. 1985; Cardelli, et al. 1991a,b; Meyer, et al. 1994, 
1997, 1998), so they can serve as reasonably good indicators for Hi in a partially ionized 
region. Sulphur is another element that has little or no depletion onto dust grains, but its 
second ionization potential is very high (23 eV). If the photoionization rate arising from 
photons with E > 23 eV is not very large, we expect that S n should be a good indicator 
of the total hydrogen in the line of sight, both in the neutral and ionized forms. In essence, 
for N, O and S we can make use of our general understanding of their depletions to arrive 
at conclusions on the relative ionization of hydrogen. 

With our ability to identify how much of the absorptions by N I and O I can be assigned 
to Components C and D, we can differentiate between the relative fractions of H in neutral 
form in each case, after making comparisons with Dupin & Gry's (1998) values for Sil in 
these same components. In doing this, we invoke two assumptions about the depletions: 

(1) the depletions of N and O onto dust grains are the same for Components C and D, and 

(2) S is undepleted. Assumption (1) seems reasonable in the light of evidence presented by 
Meyer, et al. (1997, 1998) and (2) seems justified from the evidence summarized by Savage 
& Sembach (1996) and Fitzpatrick & Spitzer (1997). If either of these assumptions is not 
quite correct, small errors in our conclusions may arise. However, they are probably not 
much worse than the uncertainty in identifying the relative contributions of Components 
C and D in the low resolution spectra of the S II features. 

For nitrogen, we obtain a general depletion factor <5(N) = 0.55 through Eq. [l] by 
taking Af(Ni)c+D, dividing it by Af(Hi) = 2. 10 18 cm~ 2 , and then multiplying the result 
by the cosmic ratio of H to N. The same analysis for O leads to 5(0) — 0.58. These values, 
2 Throughout this paper, we use the cosmic abundances of Anders & Grevesse (1989) 




(1) 
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incidentally, are not much different from those shown by Meyer, et al. (1997, 1998). We 
now estimate the hydrogen neutral fractions in each component from the expression (for 
the case of N i) 

„(H,)MH«„ M ,= (|)_ k5 ^ (2, 

and likewise for O I. For both Ni and Oi, we obtain n(Hi)/rt(H to tai) = 0.25 for Compo- 
nent C and 0.035 for Component D. 

The high column density of Si in identified with Component D by Dupin & Gry 
(1998) presents a special problem (see Table ||). The rate constant for charge exchange 
between doubly ionized Si and neutral H is about 3.0 10 _9 cm 3 s _1 at the temperatures 



of interest to us (Gargaud, et al. 1982). From the electron density derived later in §4.3 
and n(Hi)/n(H to tai) for Component D, we find that n(Hi) ~ 0.005 cm -3 (for an appli- 



cable temperature that we derive in § 4^4 ) . This density times the charge exchange rate 
constant is considerably larger than the expected ionization rate r(Sin) = 5. 10 _14 s _1 
that we infer from the hydrogen ionization and the shapes of the EUV spectra of e and 
P CMa (Vallerga & Welsh 1995; Aufdenberg, et al. 1999) which should dominate the 
photoionizing radiation field. It is hard to reconcile this inequality with Dupin & Gry's 
finding that A^(Si in) /A^(Si n) > 2 for Component D. We propose that this contradiction 
can be overcome by having the Si in exclusively within a different part of the region, 
perhaps one that is much closer to some ionizing source or, alternatively, within the 
cloud's conduction/evaporation front at an interface with gaseous material at a much 
higher temperature. There is empirical evidence that the association of Si in with lower 
ionization stages is not unusual. For instance, Cowie, et al. (1979) found that the velocity 
endpoints of Ovi, Si in, and Nil features were mutually correlated, and they used this 
evidence to suggest that the interfaces between cool and hot gas were conspicuous in the 
interstellar absorption lines. 

In sections that follow, we will ignore the existence of the Si in-bearing region and 
propose that the contradiction that we have noted justifies our regarding it as unrelated 
to the gas that holds most of the lower ionization states. 

4.3. Electron densities from N (Cu*)/N (Cu) 

The relative populations of the fine-structure states of C II are governed by the balance 
between collisions and the radiative decay of the upper level. If electrons are the dominant 
projectiles for excitation and de-excitation, the rate coefficient for de-excitations is 
8.63 10- 6 »i,2 3 _! 

72a = ^T^ C (3) 

(Spitzer 1978, p. 73), with the reverse rate given by detailed balancing, 7^2 — 
(92/91) exp(— Ei t 2/kT)j2,i- The statistical weights of the levels are <?i = 2 and g 2 = 4, 
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and the temperature equivalent for the difference in energy levels i?i )2 /fc = 
we find that the condition for equilibrium, 



94.9K. Thus 



ra(e)7i,2«(C n) = [n(e)72,i + A 2 ,i]n(C II*) (4) 
will lead to an equation for the electron density 



n(e) 





"n(ClI*) 
n(ClI) 




8.63 10- 6 fii, 2 {(f) 


cxp( fc ^ 2 ) 




'n(ClI*)~ 
. "(ClI) _ 


} 



(5) 



where the radiative decay probability for the upper level is A 2 .i = 2.29 10~ 6 s~ 1 (Nuss- 
baumer & Storey 1981), and the collision strength f2i :2 = 2.81 (Hayes & Nussbaumer 
1984). Optical pumping of the Cn fine structure levels is unlikely to happen under the 
circumstances where there is a very high optical depth in the line (Sarazin, et al. 1979). 
Also, the density of pumping radiation must be very large - even larger than the ele- 



vated levels considered later in §4.5. Finally, unacceptably high values for n(Hi) would 
be needed for collisions by neutrals to have any importance (Keenan, et al. 1986). 

As is clear from Fig. [j], the Cn feature at 1036 A is far too saturated to allow a 
determination of TV(Cli) for either Components C or D. Thus, our determination of the 
column densities must be indirect. A good surrogate for carbon is sulphur. Our repeat 
of calculations of the type done by Sofia & Jenkins (1998) indicate that in partially 
ionized gases these two elements have about the same fraction of atoms elevated to 
higher (unseen) stages of ionization for a wide range of conditions. The depletion of 
carbon atoms in dense clouds is typically 5(G) = 0.39 (Cardelli, et al. 1993, 1996; Sofia, 
et al. 1998), and, learning from the example of r CMa (Sofia, et al. 1997), this level 
of depletion seems to hold also for low density lines of sight. On the assumption that 
carbon toward (3 CMa is depleted by this factor and there is no depletion of sulphur, 
we can arrive at iV(Cn) from the product (C/S) COS mic^(C)iV(S n). Doing so gives us the 
values iV(Cn) = 8.3 10 14 cm~ 2 for Component C and 1.7 10 15 cm~ 2 for Component D, 
leading to iV(Cll*)/iV(Cll) = 0.029±£;^ and 0.0071±g;ggl| for Components C and D, 
respectively. 

Before we can derive values for n{e) from Eq. ^, we must determine T. To do this, 
we rely on another method for measuring n(e), but one that has a different temperature 
dependence. The balance between the ionization of Mg I and the recombination of Mg II is 
a fundamentally different process from that which governs the fine-structure equilibrium 



of Cn, but both are driven by the value for n(e). In the next section (§4.4), we shall 



make use of the difference to constrain other free parameter, the temperature T. 
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4-4- Electron densities from Mgi/Mgii 

The equation for the equilibrium of the lowest 2 ionization levels of Mg is given by 

[r(Mgl) + C(Mgn)n(iJ+)]n(Mgi) = 

a(Mgi)n(e)n(Mgll) (6) 

For the charge exchange rate C(Mgn) that applies to the reaction Mgi + H + — » 
Mgn + H, we used the analytical approximation C(Mgn) = 1.74 10~ 9 exp(-2.21 10 4 /T) 
derived by Allan, et al. (1988). Collisional ionization of Mgi is negligible at the temper- 
atures of interest in this study. We computed r(Mgi) at the position of the Sun using 
mean of the stellar radiation fields estimated by Jura (1974) \F\ = 3.0 10 _3 erg cm _2 s~ 1 
and Mathis, et al. (1983) XF\ = 2.5 10 _3 erg cm _2 s , multiplied by the photoioniza- 
tion cross section given by Verner, et al. (1996), to arrive at T(Mgi) =6.1 10 _11 s _1 . 
(However, this number increases to 9.5 10 _11 s _1 when we consider locations closer to 
e and (3 CMa at a later stage of the analysis.) For a(Mgl) we used the radiative and 
dielectronic recombination rates given by Shull & van Steenberg (1982), supplemented by 
the additional contributions from low-lying resonance states computed by Nussbaumer 
& Storey (1986). 

It is immediately evident that the Echelle B observations of Mg I and Mg n by Dupin 
& Gry (1998) both show an excellent signal-to-noise ratio, a smooth stellar continuum 
flux, and absorption features that are reasonably well resolved. Nevertheless, there is a 
difference between the two. On the one hand, the 2853 A feature of Mgi represents an 
absorption that is unsaturated and thus straightforward to analyze. On the other hand, 
the two features of Mgn (2796 and 2803 A) were heavily saturated, which leads to a 
large uncertainty in iV(Mgn). Even though Dupin & Gry could use the well-determined 
b value for Mgi to help in the interpretation of the Mgn features, the error in iV(Mgn) is 
still significant (see Table ||). We will use here a less direct method of deriving iV(Mgn) 
from similar but more accurately determined species, along with an application of some 
empirical relationships for interstellar depletions, that will ultimately lead to a completely 
independent estimate for iV(Mgii). 

As with C and S discussed in § |4.3| , Mg and Si are a good pair of elements that have 
very similar photoionization properties (although the rate coefficients for charge exchange 
of their doubly ionized forms with neutral hydrogen are very different - see below). 
Fitzpatrick (1997) shows a general relationship between the logarithms of the interstellar 
depletions of Mg and Si. If we assume that S is generally undepleted, we can substitute it 
for H and then derive the quantity log £(Mg) from our observed log <5(Si) and his best fit to 
the trend log<5(Mg) = 0.82 log £(Si) -0.17 dex. When we do this, we find that logd(Mg) = 
—0.63 dex for Component C and —0.80 dex for Component D. Products of these derived 
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6's, the Sll column densities, and the cosmic values for Mg relative to S give iV(Mgn) = 
5.4 10 13 cm -2 and 7.6 10 13 cm -2 for Components C and D, respectively. A moderate 
adjustment in iV(Mgii) for Component D must be made to account for the fact that 
about 24% of the Mg is probably in the doubly ionized form, while the fractional amount 
of Si II is much lower because it has a significantly larger rate coefficient for recombining 
via charge exchange with neutral hydrogen. This conclusion was derived quantitatively 
after making preliminary calculations of the ionization balance, as described later in 
§ [i~5| . No such correction is needed for Component C. The estimate we made here for 
iV(Mgn) for Component C is somewhat larger than the one given in Table || - however 
is marginally consistent with it when the error of 0.16 dex attached to our present method 
is considered- while our (adjusted) value of iV(Mgll) = 5.8 10 13 cm -2 is remarkably close 
to the value derived by Dupin & Gry. Finally, we arrive at iV(Mgi)/iV(Mgii) equal to 
O.OQ15±g;|$8£ for Component C and 0.0043±£^ for Component D. Note that while 
Dupin & Gry's values for iV(Sn) have large estimated relative errors (36%, 22%), such 
errors have a very weak influence in the conclusions because the slope of the empirical 
relationship between logJ(Mg) and log (5 (Si) is not much different than 1.0. The same 
applies to errors in the assumption that log<5(S) = 0.0 dex. 

As a check that the depletions toward (3 CMa are not out of the ordinary, we carried 
out a similar exercise, comparing the measured log<5(Fe) = (—1.20, — 1.60) dex and 
their respective logJ(Si) = (—0.56, — 0.77) dex for consistency with the trend shown by 
Fitzpatrick (1996). Both cases fell within the body of points that defined this trend. 

The curved bands in Fig. || show, as a function of temperature T, the expected out- 
comes for log AT(Mgi) — logAT(Mgn) found by solving the ionization equilibrium equa- 
tion (Eq. ||) with the values of n(e) obtained from Eq. ^| and the measured values of 
AT(Cll*)/iV(Cil) given in § }4.3| . From the intersections of these curves with the loga- 
rithms of the observed iV(Mgi)/iV(Mgii) presented above (straight, horizontal bands in 
the figure), we find that 400 < T < 6500 K for Component C and 8000 < T < 14 000 K 
for Component D. To compute the upper limits for n(e), we evaluate Eq. || with the upper 
limit for N (C II*) / TV '(C II) at the higest temperature that is consistent with this upper 
limit (not the highest temperature allowed in general). For the lower limits, we take the 
lowest permissible iV(Cll*)/iV(Cii) and apply it to Eq. || with the lowest temperature. 
These conditions for the limits are shown by solid dots in Fig. ^. For Component C we find 
that 0.08 < n(e) < 0.6 cm~ 3 , and for Component D we arrive at 0.09 < n(e) < 0.2 cm~ 3 . 
(While our lower limit for n(e) for Component C is formally allowed, the real value of 
n(e) is probably much closer to the upper limit because the temperature of the gas is 
probably much higher than 400 K.) These represent the worst possible extremes in n(e) 
permitted by the data. 
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Fig. 2. Predictions for the logarithms of 7V(Mgi)/iV(Mgii) as a function of temperature 
T (curved, shaded bands), computed from Eq. |g| using electron densities n(e) derived from 
the measurements of iV(Cll*)/iV(Cll) applied to Eq. 0. The widths of these bands indi- 
cate the ranges that are allowed by our estimated errors. The shaded, straight, horizontal 
bands indicate the observed values (and their uncertainties) for log[iV(Mg i) /iV(Mg n)] . In 
each case, the overlaps of the straight and curved bands indicate the temperature ranges 
that are consistent with both the N(C II*)/JV(C II) and iV(Mgi)/iV(Mgii) results. Po- 
sitions on the diagram that represent the worst extremes in n{e) are indicated by solid 
dots; sec the text for details. 



4-5. Ionizing radiation 

Shortward of the Lyman limit, the general ionizing radiation field near the Sun is domi- 
nated by the output of e and (3 CMa supplemented by many nearby white dwarf stars, 
cataclysmic variables, and the line emission from some late-type stars (Vallerga 1997). 
In addition to stellar sources, some diffuse energetic radiation from the hot gas that sur- 
rounds the LIC may contribute additional ionization (Cheng & Bruhweiler 1990). The 
intensity of this diffuse radiation is rather uncertain: Jelinsky, et al. (1995) have observa- 
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tions that place a meaningful upper limit for this radiation if it produces spectral features 
in agreement with some specific predictions of hot-gas emission models. 

If we calculate the hydrogen ionization produced by Vallerga's (1997) composite EUV 
stellar radiation field, supplemented by Cheng & Bruhweiler's (1990) hot gas radiation 
attenuated by absorption from a hydrogen column iV(Hi) = 2 10 17 cm~ 2 , we find that 
the neutral fraction of hydrogen is 0.89 and 0.47, for Components C and D, respectively, 
if we require that n(e) — 0.31 cm -3 and 0.13cm~ 3 as we found from the analysis of 



7V(Cii*)/iV(Cli) and iV(Mgi)/iV(Mgii) in §§ 4.3 and 4.4. This is clearly inconsistent 



with our finding iV(H i)/7V(H to tai) = 0.25 and 0.035 for the two components in §4.2| . The 
model calculations followed the scheme outlined by Sofia & Jenkins (1998) that included 
the ionization of helium and allowed for the extra electrons coming from He. (In contrast 
to the LIC, where helium ionization is important, the models indicated that the high 
values of n(e) effectively suppressed the ionization of helium.) 

To overcome the fact that the predicted ionization of H is less than observed, we must 
find additional sources of ionization. One possibility is that there is vestigial ionization 
left over from a previous event that heated the gas and collisionally ionized it (Lyu & 
Bruhweiler 1996), such as the passage of a shock front in less a few times 10 5 yr ago. If this 
were the case, it would be difficult to explain why the gas is not moving at a substantial 
velocity. Another tactic is to propose that most of the gas that we see is rather close to 
e and CMa, where the hydrogen (and helium) ionization rates must be much higher 
than in the local vicinity. 

We experimented with ionization models that allowed for the enhancement of r(Hi) 
by placing the clouds closer to CMa. For one possible location, at the nearest point to 
a CMa with its white dwarf companion, the increase in the radiation field is not much: 
the strength of the radiation from Sirius B (Holberg, et al. 1998) reaches a maximum of 
only 0.4 times that of CMa as seen near the Sun. However this applies only to fluxes 
observed at wavelengths above the Lyman limit. If there were substantially less Hi and 
He I between Sirius B and the cloud in front of CMa, the ionizing fluxes could be larger. 

We expect that at a distance of 130 pc from the Sun, a cloud is at the smallest 
possible distance from e CMa (33 pc) and only 26 pc from CMa (assuming that the 
most probable distances between us and the stars derived from Hipparcos are correct). 
Under this condition, the flux of e CMa is enhanced by a factor of 17, and the flux from 
CMa increases by a factor of 36. With the spectral distributions given by Vallerga 
k Welsh (1995) for e CMa and Aufdenberg, et al. (1999) for CMa (after correction 
for absorption by the H and He between these stars and the Sun), we find that the 
enhanced fluxes at this position produce hydrogen neutral fractions of 0.25 and 0.037 for 



Components C and D. These values are close to what we observed and reported in §4.2. 
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To some degree, our results are dependent on the assumed iV(Hi) and iV(Hel) be- 
tween the ionizing sources and most of the gas in the components. Absorption in front of 
(3 CMa is constrained by the column densities that we observed, whereas those for e CMa 
are arbitrary. For Component C our adopted columns were iV(Hi) = 1.1 10 18 cm~ 2 for 
(3 CMa and 1.1 10 17 cm~ 2 for e CMa (this low value was required to raise the ionization 
high enough). For Component D, we made the columns in front of both e CMa and 
j3 CMa equal to 7V(Hi) = 2 10 17 cm -2 . For both components, we made the Hei to Hi 
ratio in the absorbing gas consistent with what we found from the ionization equilibrium 
calculations. 

Our ionization equilibrium calculations indicate that both clouds are located in a 
region close to e and f3 CMa where the radiation fields are strongly enhanced. Note 
also that a physically separate region, perhaps an outer layer of the cloud representing 
Component D that has evolved to a more fully ionized state, is revealed by the presence 
of Sim. 

The fractional ionization of He that we compute is not very substantial. Singly ionized 
He accounts for only 6% (Component C) and 19% (Component D) of the total He, even 
when the material is placed at the position where the radiation fields from e and (3 CMa 
are enhanced by the large factors given above. The total column density iV(Hei)c+D = 
1.6 10 18 cm~ 2 , which is above the lower limits stated in jj |4.2| . For Component C, the 
computed Hei/Hi of 0.38 agrees with the general results of Wolff, et al. (1999), while the 
much higher ionization of H in Component D produces a considerably larger He i/H 1=2.2. 
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